Quantum Spin Dynamics of Spin-1 Bose Gas 
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We show that the quantum evolution of a spin-1 Bose gas with nearly all bosons initially in 
the F z = state has a "quantum carpet" spin-time structure with self-similar properties. The 
system continuously evolves into "multi-peaked" Schrodinger cat like states, returning occasionally 
to coherent structures, which leads to large number fluctuations (as seen in recent experiments). 
The self similar behavior allows one to reveal the quantum evolution as a set of peaks in the number 
probability distribution of a spin component at times much shorter than the quantum revival time. 
We also show that these features survive small number fluctuations among spin components up to 
a few percent. 
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Since the discovery of Bosc-Einstein condensation in 
alkali atomic gases, there have been increasing activities 
in Bose gases with non-zero spin F. Compared to scalar 
Bose gases, these systems have a greater variety of ground 
states and macroscopic quantum phenomena. Bosons 
with spin also show great promise for quantum compu- 
tation and quantum information processing, since their 
spin-dependent Hamiltonians have been derived and are 
completely specified by the scattering lengths in differ- 
ent angular momentum channels Q. Despite the sim- 
ple appearance of these Hamiltonians, their dynamics is 
surprisingly rich. Understanding spin dynamics, and in 
particular quantum spin dynamics, is important not only 
for theoretical reasons, but also for engineering and con- 
trolling the time evolution of spin carrying bosons. 

Recently, the spin dynamics of F = 1 and F = 2 
Bose gases of 87 Rb have been experimentally studied . 
For F = 1 87 Rb bosons, if initially all bosons are in 
the F z = level, the number of bosons in this level at 
later times (when averaged over many runs) exhibits a 
"damped" oscillatory behavior with very large fluctua- 
tions, whereas the fluctuations in the number difference 
between F z = 1 and F z = — 1 spin states is almost zero 
(due to spin conservation). At present, the theoretical ex- 
planation for this oscillatory behavior is based on mean 
field spin dynamics within the so-called single mode ap- 
proximation (SMA), with an averaging over the initial 
condition to simulate the slightly different initial condi- 
tions in different experimental runs. (SMA means all spin 
states have the same spatial density profile). The large 
fluctuations are due to the sensitive dependence of the 
semiclassical evolution on the initial conditions Q . Mean 
field spin dynamics assumes that the system remains an 
un-fragmented condensate at all times. This assumption 
is considered satisfactory because the mean field calcula- 
tions seem to account for the observed damped oscillation 
in the F z = component. There is, however, no direct 
experimental verification of it. 

On the other hand, it is known that quantum dynamics 
also exhibits a "damped" oscillation with large fluctua- 
tions in the F z = component 0, although this oscilla- 
tion is in fact non-dissipative with a period t a . In current 



experiments this quantum revival time is longer than the 
duration of the experiments and is thus not observable. 
In this letter we present a method for distinguishing be- 
tween the semiclassical and quantum dynamics in time 
scales much shorter than t a . We shall show that hidden 
in the quantum spin dynamics of spin-1 Bose gases lies 
a set of beautiful scaling relations in the noise statistics. 
Not only are they dramatic demonstrations of quantum 
behavior, but also practical diagnostics that allow us to 
detect the existence of quantum revivals at times smaller 
than the full revival period. 

(A) Mean field spin dynamics: Within SMA, the 
Hamiltonian of N spin-1 bosons in zero magnetic field 
isBUl H = cS 2 , S = Y^pv ojiS^Oy, where creates a 
boson in the state F z = ji and Su V is the spin-1 matrix. 
We can also write S 2 = {N x - N_i) 2 + (2iV - l)(Ni + 
N-i)+(2a\a_ 1 a 2 ,+h..c.), where N M — aj^a^ is the number 
of bosons in state \i. The Hilbert space is spanned by 
the states \N X , N , N_ x ) which have bosons in state 
fi. In the presence of a magnetic field B = Bi, there 
are linear and quadratic Zeeman effects given by H± = 
-p(Nx - JV_i) and H 2 = a(N$ + iV 2 ^), where p and a 
are proportional to B and B 2 respectively. Due to the 
conservation of the total spin S z = M, the dynamics is 
independent of p. Since we are interested in low field 
phenomena, we will consider the case a = 0. 

The mean field approach assumes that the system is in 
a coherent state |*(t)) = (£ M Cm (<)oJ,) JV |0}/\/M, where 
Cfj.(t) is a normalized vector. The single particle density 
matrix is = (*(t)|aj, a„|¥(t)) = NQ(t)CJt), which 
has only one eigenvalue. The occupation number of state 
H is N^t) = |$ M (t)| 2 , where = y/N^ft) is the 

condensate wavefunction which evolves according to the 
mean field equation ifid^^t) = 2c(S(i)) • S„„<& !/ (t), and 
(S(*)) - Eo/j^WS^^W. If initially all the bosons 
are in the fi = state as in ref. 0, (i.e. 3^(0) = VNS^fi) 
the system does not evolve, since the mean field equations 
imply that such a state is stationary. In practice it is 
assumed that in each run of the experiment there is a 
small but non-zero $±1(0). The experimental results are 
to be compared to the average of such runs, which yields a 
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FIG. 1: Probability of finding p bosons in F z — +1 state as a function of time r for 1000 bosons. Light (dark) colors are 
regions of large (small) probability. At r = (r = 1/2), the probability is sharply peaked at p = (p — n). In between, the 
distribution of probabilities has a fractal behavior. 



decay of the oscillations in No(t) with large fluctuations. 

(B) Quantum spin dynamics: In contrast, in the 
quantum dynamics the state |0, iV, 0) is not stationary, 
due to the "conversion" term aja^a 2 , in S 2 . In the 
following, we shall focus on the evolution of this state, 
since it is related to the experimental situation, and is 
the simplest case that captures most of the intricacy and 
richness of the quantum evolution. For simplicity, we 
shall consider a system with an even number of bosons, 
N = 2n, where n is an integer. Since H a is rotation- 
ally invariant, it only couples the initial state |0, N, 0) 
to states \p) = \p, N — 2p, p) , where p is the number 
of bosons in state /i = ±l, (0 < p < n). We can 
then expand the quantum state at time t as |^(t)) = 
£p=0*(P.*)b>> = (ple-u^lO), with initial 

condition 0) = 5 Pt Q. For operators O that are di- 
agonal in number space, their expectation values are 
(0(t)) = E p O(p)\^(p,t)\ 2 . For N = 2n » 1, it is 
convenient to write 



p/n = (cos 6 + l)/2, < 9 < vr, 



(1) 



*(p,t) = tf(0,t), s o that (6 (t)) = ^d9 0(9)\f(9,t)\ 2 , 
where f(9,t) = ^(n/2)sin6» *(6», t). The key features 
of the quantum dynamics is contained in the behav- 
ior of f(6). The system will be either a number (i.e. 
Fock) state or a coherent state if the probability distri- 
bution V{9, t) = \ f(9, t)\ 2 is a delta-function or a sharply 
peaked function of width y/n, with a number fluctua- 
tion (AiV ) 2 = ((N - (N )) 2 ) ranging from to n. If, 
however, f(9,t) is a sum of several sharply peaked dis- 
tributions, the system is in a Schrodinger cat like state, 
with (AN ) 2 of order n 2 . 

(Bl) Quantum Carpet: Since H a = cS 2 , 



*(p,t) = J2(p\s,o)(s,o\o)( 



(2) 



where r = t/t a , t Q — nh/c (independent of number of 
particles), \S, M) is the state with total angular momen- 



tum S, and S z — M. Since N is even, Bose statistics 
implies that S is also even, and we shall write S = 2k, 
where < k < n. Since (p\S, 0) is real (see Appendix), 
it is clear from eq.J2J 

/(0,t)=/(0,t + 1), f(0,-T) = f(6,T)*. (3) 

Using the exact expressions for (p\2k, 0) derived in the 
Appendix, we can calculate the probability | ^ (js, t) | 2 ex- 
actly for any number of particles. In figure 1, we display 
\^(p, t)| 2 in a "spin-time" plot for a system of 1000 par- 
ticles, where light color means high probability. It is a 
highly intricate yet orderly self similar pattern, analo- 
gous to the so-called "quantum carpet" that displays the 
motion of a quantum mechanical particle in a one dimen- 
sional box@. At r = 0,1/2, and 1/4, |*(p, r)| 2 consists 
of a single sharp peak at p = 0, n, and n/2 respectively. 
They are therefore either a Fock state or a coherent state. 
At any other time, however, \^f(p, t)| 2 has many different 
peaks and the system is therefore Schrodinger cat-like. 
For example, at time r = 1/8 and 1/3 one can see from 
figfflthat the distribution consists of two sharp peaks. 

To understand these intricate features we need to ana- 
lyze the spinor wavefunction and their initial probability 
amplitudes = (2k, 0|0). In the appendix we show that 
for 1 <C k <C n,p, (n — p), the spinor wavefunction can be 
approximated by (up to 0(l/n)) 0, 



(p|2fc,0) 



4 (-l)P+ k 



cos 



(k + \)9 



Ttn Vsin6> 

and that to the same accuracy, we also have 
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4k + 1 
2n 



-k 2 /(2n) _ Ae (k~k f/o 



; (4) 



(•5) 



where k ~ 0. 8y/n, a ~ y/n, and A is a constant. Thus 
the initial state populates angular momentum states up 
to S of order y/n and we are allowed to use We then 
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have 



f(9, r)=J2\ ±Ck e - l2 ^ 2k+1 ^ +l ^P+ k ) cos 



{k + - A )0 



(6) 

It proves useful to extend / outside the physical inter- 
val [0, tt] to all 9, so that f(9, r) = f(9+8n, r). Eq.© im- 
plies that f(9, 0) is a sequence of sharp peaks at (21 + l)ir 
(where I is an integer) such that f(iv,0) — —f(3n,0) — 
—f(5n,0) — +f(7n,0). It is straightforward to show 
from eq.© that 



f(0,T 



E 



f{0,T+ 1 -) = (/(0-7T,T)+/(0 + 7r,T))/V2. (7) 

Finally, eq.© also implies (shown in next section) that 
if 9/(4t + 1) < 07 



JC{r) 



4r + 1 



4r + 1 ' At + 1 



(8) 



where C(t) = 3 ^ t (t 2 + (0/2^) 2 ). 

We are now ready to study fig.l. Consider the dense set 
of rational times r = P/Q, where P and Q are integers 
relatively prime to each other, and note that 

(a) : Eq.© implies any rational r = P/Q > 1/2 can be 
"mapped" onto P'/Q < 1/2 with P' < P. It is thus 
sufficient to consider rationals in [0, 1/2]. 

(b) : Eqs. © and © imply P/Q e (0,1/2) (except for 
1/4) can be mapped onto r' = P/Q' with Q' = (4P-Q) 
such that |Q'| < Q. 

(c) : Since /(0,O) is sharply peaked at (2£ + 1)tt, Eq.© 
implies f (9,1/2) peaks at 2^7r and /(#, 1/4) peaks at (£ + 
l/2)7r. The latter is due to the fact that f (9,1/4) is a 
sum of /(#, 0)'s moved forward and backward by 7r/2. 

Since both (a) and (b) reduce P and Q, their repeated 
application can reduce the distribution at r = P/Q to 
that at Tp = 1/Qf = 1,1/2, or 1/4, depending on 
whether Q is odd, of the form 4m — 2 or Am respec- 
tively. Using the fact that after one application of eq. © 
the 9 variable gets scaled by a factor Q'/Q, successive 
application of (a) and (b) to the final stage will scale 9 
by a factor Qf/Q- This leads to the following "Rules" 
for the peaks 9* of the distribution function. 

(1) if Q = 2m- 1, Q F /Q = 1/Q, then 9* = k(2£+1)/Q, 

(2) if Q = 4m - 2, Qf/Q = 2/Q, then 9* = A£ir/Q, 

(3) if Q = Am, Q F /Q = 4/Q, then 9* = (2 + A£)ir/Q; 
where £ = 0, ...,m— 1. In each case there are m peaks in 
the physical interval [0, 7r]. 

As an example of the application of these rules, we list 
the peak locations of the probability \^(p, P/Q)\ 2 : 
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Case 
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9* 


A6»* 


p* /n 
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(2) 
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1/2 
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(1) 
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7T/0, 7T 


27r/3 


1/4, 
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(3) 


1 


tt/2 




1/2 


5 


(1) 


3 


7t/5, 37T/5, 7T 


2vr/5 


eq.JTJ 


20 


(3) 


5 


for/lO, fe = 1,3,5,7,9 


tt/5 


eq.lJTJl 



The case r = 1/20 is shown (as the black curve) in the 
main part of fig. |3 

The intricate behavior in figUis due to the fact that 
as r increases, it sweeps through a dense set of rationals 
where Q changes in a chaotic fashion. The important 
point, however, is that even when the quantum revival 
time is longer than the duration of experiment, its pres- 
ence is reflected as a sequence of distinct peaks in the 
probability distribution at shorter times. The observa- 
tion of these peaks in the probability distribution would 
be a clear demonstration of quantum evolution. 

(C) Derivation of eq.©: Recall that the elliptic 
theta-function 9 3 (u,t) = J2 k (e l7TT ) k2 e 2i7Tkv obeys 



Q 3 (v,t) = B 3 (v,t + 2) =G(v 
and satisfies the reciprocal relation^ 



Q 3 (v,T) = J-e-^ v ' r Q 



1,70 



(9) 



(10) 



It is then straightforward to show that eq.© is 

f(9,r)= B(-l)^ ± e W/4 

x9 3 (±£ -T- ik e,-(Ar + 1) + ie) (11) 

where e = l/(7rcr 2 ) and B is a constant . To obtain the 
scaling relation © we apply (|10fl to both theta functions 
in ea. Hll(l . neglecting the imaginary parts in the expo- 
nential parts. This is valid as long as 9k e -C (At + 1). If 
we further use © to subtract 2 from the second variable 
of the Theta function and use that Q(—v,t) — Q(v,t), 
we obtain ©. The restriction on the imaginary parts 
shows that the maximum value of Q for which the rules 
can be applied is of the order y/n. 

(D) Fluctuations in initial states: Since the semi- 
classical evolution of states close to |0,iV, 0) is highly 
sensitive to initial condition, we need to consider how 
the quantum evolution is affected by small fluctuations 
in particle numbers in various spin components. The 
state of interest is now \^(t)) = J2 P ^(Pi *)bi -^0, where 
\p; M) = \p + M, N - M - 2p, p) and M is the mag- 
netization. The initial condition is |^(0)) = \q\M) 
with q,M « N, and we are interested in \&(p, r) = 
J2 k (p; M\S, M)(S, M\q; M)e- iTrS ^ s+1 ^ . We present here 
the numerical results and describe t he p henomena but we 
shall present the details elsewhere [llj. 

Let us consider the case M = Q, m = n_i = q. When 
q =/= 0, each of the m peaks in the original probability 
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FIG. 2: Probability of rinding p bosons in the F z — 1 state at 
t = 1/20 for N = 2n — 1000. The black, red, and green lines 
correspond to the initial state \q, 2(n — g), (?) with g = 0, 4, 8 
resp. The inset shows the probabilities for r = 1/4 for g = 
0,4,8,20,40 (black, red, green, blue and black lines, resp.). 



distribution \f(9, P/Q)\ 2 (separated by A9* in 0-space) 
is split into a group of q + l-subpeaks centered around the 
original peak at 9*, such that the outermost two peaks 
have the highest probability. The range of the splitting, 
i.e. the distance between the two outermost peaks, is 
59 q = Ay/ q/n 11]. In order for the original peaks to 
remain as separate groups, one needs 89 q < A9* . The 
critical fluctuation (q/n) c beyond which different groups 
begin to interfere is (q/n) c = (A9*/4) 2 . 

This phenomenon is demonstrated in fig0 where we 
have plotted \^(p, r)| 2 for a system of 1000 bosons with 
t = 1/20 for different number fluctuations q = 0,4,8, 
corresponding to q/n = 0, 0.008, 0.016 respectively. In 
this case Q = 20, A9* — ir/5, and we have (q/n) c — 
0.025. For q = 4, q/n < (q/n) c , we see that while each 
original peak is split into five subpeaks, each group of 
peaks is still separated from the rest. Another example 
of the effect of q on the original peaks is shown in the in- 
set for the probabilities at r = 1/4 for q = 0, 4, 8, 20, and 
40, with q/n ranging from to 8%. Given that differ- 
ent values of q are realized in different runs, experiments 
sample the average of these curves. What is remarkable 
is that for the case of r = 1/4, the original single peak 
structure in the probability distribution is still reflected 
in the averaging process for q/n as high as 8%. These 
behavior also applies to much larger systems, as well as 
states with M^0. 

We have therefore demonstrated that the peak struc- 
ture in the probability distribution at rational times, 
which is the hallmark of quantum evolution, survives 
number fluctuations in different spin components as long 
as they are no more than a few percent. 

Appendix: The fully polarized spin states \S = 
2k, 2k) can be expressed as 12] a\ 2k e t "- fe |0)/ v //(n; k) 
where the normalization constant is /(n; k) = (2fc)! 2 (n — 
k)\ (2ra + 2fc+l)!/[(n+fc)!(4fc+l)!], and 9 = a§-2aia_i 



is the singlet pair destruction operator. The states 
with the same spin but M = are obtained by ap- 
plying the S_ operator 2k times, yielding \2k, 0) = 
5 ,2fe |2fc,2fc)/vT4fc)! with S_ = v / 2(aJ°i + a -i a o)- Since 
[0,5-] = 0, we can permute the order of these 
two operators, and we are lead to consider (ajai + 



aLjao) 



2fc a t2fc, 



= Ei=o(( a o 2fe Ja o)) a "i a i 



ti 2k-jJ2k 



a 



where we define {{b l c 3 )) as the sum of all different or- 
ders in which the i + j operators can be arranged. 
We can show by induction that ((aj 2 _ ^ c 



mi 



(2k-2j)\2ij\ u o 



t2(fc-j) 



|0) and binomially expanding 8^™ 
(p|2M> = A(k,n)p\^(2n-2p)\^ZZ n(k ' n - p) 

(-4) 3 (2fc)! 2 (ra-fc)! 
(2s)! (k-s)Y 2 (n-p-s)! (p-fc+s)! ' 

approximation we get (JSJ) . 



max(0,fc- p) 

For p — using Stirling's 



To verify eq. 
N\N S for s < 



we use that (N + s)\ w 
N, and thus obtain (p|2fc,0) w 

(4fc + l)(2n-2p)T (-2)P(2fc)! 2 «! j? ,_, l.i _ H \ 

2n(2n)! (-4) fc (n-p)! k\ 2 V n ) 2 1 V "'J ^f^' p> 

where ofi is the hypergeometric function. Using now 



that |9J 2 F l ia,(3- T ,z) = (1 - z)- a 2 F 1 (a, 7 - 0; 7; ^) 



and that 



FUa, 
M 



iFi(-k,k + 1 + a + 0;a + l;x) = 



r[l { + T+k) P^i 1 -^) wh ere P[ a ^\z) is a Jacobi poly- 
nomial, and finally using that for k 3> 1 (although still 
k « n, P ) Pi-^W) - v , fc : os( , /2) cos [(* + 
as well as the Stirling approximation we arrive at (|4"|l. 
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